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Abstract. We give an outline of the computation of the third order term in 
a generaUzation of the Strong Szego Limit Theorem for a zeroth order pseudo- 
differential operator (PsDO) on a ZoU manifold of an arbitrary dimension, see 
[ |Gi2| for the detailed proof. This is a refinement of a result by V . Guillemin and 
K. Okikiolu who have computed the second order term in [ G102 1 . An important 



role in our proof is played by a certain combinatorial identity which generalizes 
the f ormula of G. A. Hunt and F. J. Dyson to an arbitrary natural power, see 
[ |Gi3[ . This identity is a different form of the renowned Bohnenblust— Spitzer 
combinatorial theorem which is related to the maximum of a random walk 
with i.i.d. steps on the real line. 

A corollary of our main result is a fourth order Szego type asymptotics 
for a zeroth order PsDO on the unit circle, which in matrix terms gives a 
fourth order asymptotic formula for the determinant of the truncated sum 
Pn{T\ + T2D)Pn of a Tooplitz matrix T\ with the product of another Toeplitz 
matrix T2 and a diagonal matrix D of the form diag(- ■ • , i , i , 1, 0, 1, ^ , ^ , ■ ■ ■ ). 
Here P„ = diag(- ■ • , 0, 1, • ■ ■ , 1, 0, ■ ■ ■ ), (2n + 1) ones. 



1. Introduction and main results 

The main motivation for this work was to find an explicit formula for a "Szego- 
regularized" determinant of a zeroth order pseudodifferential operator (PsDO) on 
a ZoU manifold introduced in [GOl, p2| , see Remark 2.4. 
is valid for any dimension d G 



1.2 



N. 



Our main result 
In the case d = 2, Theorem 1.2 



The- 



gives 



such a formula. 

In this paper we find a third order generalization of the Strong Szego Limit 
Theorem (SSLT) for a zeroth order PsDO on the unit circle (Theorem 1.1), and on 
a Zoll manifold of an arbitrary dimension d e N (Theorem ^]^). We give also an 
outline of the proof of Theorem 1.2 and derive Theorem |l.l| from the former. The 
detailed proof of Theore m |1.2| can be found in | Gil , Chapter 1] and will be published 
in a forthcoming paper |Gi2[ . In Section 3, we sketch a proof of the generalized 

for the detailed proof. 



Hunt-Dyson combinatorial formula (Theorem 3.1), see |Gi3 
We give also a brief review of the related work. 

Recall that M is called a Zoll manifold if it is compact, closed and such that 
the geodesic flow on M is simply periodic with period 27r. The unit circle and 
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the standard sphere of any dimension are ZoU manifolds. A second order gener- 
alization of the SSLT for a ZoU manifold of any dimension has been obtained by 
V. Guillemin and K. Okikiohi |G01, G02|, see also an important preceding work 
|01| by K. Okikiohi for the case of the two- and three-dimensional sphere. The 
proofs in [01, GOl, G02| use a combinatorial identity due to G. A. Hunt and 
F. J. Dyson and proceed in the spirit of the combinatorial proof of the classical 
SSLT by M. Kac |Q. See also |G03, 02 1 where the combinatorial approach and 
the Hunt-Dyson formula are used in a different setting to obtain a second order 
generalization of the SSLT for a manifold with the set of closed geodesies of measure 
zero in the unit cotangent bundle. 

In the proof of Theorem 1.2 we use the method of |G02|. A central role in 
our proof is played by a certain combinatorial identity which generalizes the Hunt- 
Dyson formula mentioned above to an arbitrary natural po wer. We call this identity 
the generalized Hunt-Dyson formula (gHD), see Theorem ^.l[ . 

After having discovered and proved the formula gHD we realized that it is related 
to another combinatorial theorem, which has a long history. The mentioned theo- 
rem is a result due to H. F. Bohnenblust that appeared in an article by F. Spitzer on 
random walks ||s|. Theorem 2.2], and is now commonly known as the Bohnenblust- 
Spitzer theorem (BSt). The characteristic function of the maximum of a random 
walk with independent identically distributed steps has been computed in ||s| with 
the help of the BSt. On the other hand, using the gHD, we can compute ||Gi3| 
the moment of an arbitrary order of the mentioned maximum (note that the usual 
Hunt-Dyson formula allows one to compute only the expected value of the maxi- 
mum, see Q). This indicates that the gHD and the BSt should be essentially the 
same. And indeed this is the case: it turns out that the gHD can be derived from 
the BSt, and vice versa, see [Gi3|. 



We have found the formula gHD being unaware of the BSt. In |Gil, Chapter 2], 
a proof of the gHD "from scratch" can be found. This together with a a derivation 
of the BSt from the the gHD in [Gil, Gi3| gives a new proof of the BSt. There 



are known at least four other proofs: the original one [|s|, a proof by G. Baxter 



Bll B2|, by J. G. Wendel |Wc[, and finally, a unifying approach of G.-C. Rota 
in the framework of (Glen) Baxter algebras [Rc, RoSm| ]. See Section 3 for more 
combinatorial details. 

Let det denote the determinant of a finite rank operator. In Section 2, we give 
explicit asymptotic formulas for logdetP„i?P„, n — > oo, where P„ is the projection 
onto the first n eigenspaces of the Laplace-Beltrami operator on a ZoU manifold, and 
B is an arbitrary zeroth order PsDO, see Corollary 2f2 (for d= 1) and Corollary 2^ 
{d> 2). 

Let S"'^ be the unit circle R/27rZ. Denote by P„, n e N, the orthogonal projection 
from L^(§^) to the subspace spanned by {e^^^}\k\<n- For a function / e L^(§^) 
denote its fcth Fourier coefficient by 



(1.1) 



fk:= 



27r 



f{x)e 



dx 
2^' 



A: G Z. 



Let h{x) be a strictly positive function on such that X^fcez 1^1 l(log^)fcP < 
Denote by B the operator of multiplication by b acting in L^(S^). Recall that the 
matrix representation of the operator B is the Toeplitz matrix (bj-k)j,kez,- The 
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classical Strong Szego Limit Theorem (SSLT) states that 

oo 

TrlogP„BP„ = TrP„(IogB)P„ + ^ A: (k^)^ + o(l), 



fe=i 



Here Tr log P„PP„ = log det P„BP„ and Tr P„ (log P)P„ = (2n + 1) " log 6( 
It has been shown by H. Widom that the remainder is 0{n^°°) if b{x) € C°°(§^), 
see 



dx 
2-K • 



Wl 



Let M be a closed manifold of dimension d e N. Let ^"^{M), m g Z, denote the 
space of the classical PsDO of order m on M . Recall that for a given G G ^'™(M), 
its principal symbol (Jm{G) and subprincipal symbol sub(G) are well-defined on 
T*M. Let S*M be the unit cotangent bundle and denote by dxd^ the standard 
measure on S* M divided by {2'kY. 

The simplest form of our result is for d — 1 and the function ^/'(u) = logu which 
is analytic in a neighborhood of u = 1. 

Theorem 1.1. Let M = and Pn he the projection on the linear span o/{e*'^^}|fc|<„. 
Let B G '^'^{M) and assume that <Jq{B) is strictly positive and a certain symbolic 
norm of I — B is sufficiently small. Then \ogB G and the following holds 

(1.2) 



Tr log P„BP„ 



TrP„(logP)P„ + - 



oo 

E 

k=l 



fc(ao(logP)),(ao(logP))_,,dxde 



1 1 
n ' 2 



S'M 



E 

k=l 



k (ao (log B))^ (sub(log P))_^. dxd^ + O 



In ( |l.2| ) the argument (x,^) G S*M is omitted for brevity. By the Fourier 
coefficient we mean the following: for a fixed (x,^) G S*M compute the Fourier 
coefficient in accordance with (IT) along the unit circle being the closed geodesic 
starting at (a;, ^). 

Let us denote bo := (Tq{B) and 6sub := sub(P). Recall that for an analytic /, 
cro(/(B)) = f{(To{B)) and sub /(P) = /'(cro(^)) sub(P). Then (O) takes the form 



(1.3) 
Tr log P„PP„ 



TrP„(logP)P„ 



> fc (log 6o ) (log 6o ) _ dxd^ 



1 1 

n ' 2 



S'M 



oo 

E 

fe=i 



k {\ogbo)f.{bsuh/bo) _^dxd^ + 0[ 



Remark 1.1. Theorem 1.1 is a particular case of Theorem 1.2 below, see Remark 1.4 



Let us fix some notations and then state the result for the higher dimensional 
case. Let M be a Zoll manifold of dimension d G N. Let A denote the Laplace- 
Bcltrami operator on M. It is known |DG] that there exists a constant a G M such 
that the spectrum of \J lies in bands around the points fc-l- fc G N. Moreover, 
in has been shown in |CdV| that there exists A-\ G 4'^^(Af) such that [A, = 
and the spectrum of the operator 



4 



DIMITRI GIOEV 



is N. Let Pn, n E N, denote the projection from L?{M) onto the subspace spanned 
by the eigenfunctions of A corresponding to the eigenvalues 1,2, •• • ,n. Let dxd£, 
be the defined above measure on S*M := {(x,^) : ai{A){x,£^) = 1}. Following 
|G01| , |G02[ , for a S G 4'"(M) and i e R introduce the operator 

(L5) := e-'^-^Be'*-^. 

By Egorov's theorem, e ^'"(M), and also 

where 9' stands for the shift by t units along geodesic flow. Note that because 
sub (A) = const the following also holds ||^, Lemma 2.2] 

snh{B'){x,0 = sub(B)(e*(a;,C)). 

Because spec(A) = N, the operator _B* is periodic in t with period 27r. 

The result we are about to state, and also Theorem 1.1 above, gives information 
on the asymptotic behavior of 

(1.6) Tri^(P„BP„)-TrP„?/.(S)P„, n ^ oo, 

for the analytic in a neighborhood of 1 function ip{u) = log u. In | Gil , |Gi2 ] we obtain 
the corresponding statements for an arbitrary analytic function 'il'{u). These results 
involve certain maps W, <&, T whose action on an analytic function is a continuous 
function from to C for some j € N. The map 



Wmyuy2) ■■= 



1 



1 ^/{u,)-i,'{u2) 



Ul — U2 



duidu2, yi,y2 > 0, 



was obtained by A. Laptev, D. Robert and Yu. Safarov in [LRS] (here j — 2). 
Earlier, an equivalent to W map has been obtained by H. Widom [ |W2| , |W3| in a 
two-term Szego type asymptotics for integral operators with discontinuous symbols. 
It has been noticed by H. Widom, and also by the authors |LRS, 02 1, that for the 
function V'(w) = logw, the action VFpog] takes a simple form 

W [log] (yi , y2 ) = - ^ log y 1 log 2/2 • 

This explains the fact that the map W is not required in the next theorem. 
Now define for j/i, ?/2, ys < 1, 



*[log](2/i,2;2,y3) = ys 



(1.7) 



yi ry2 



Jo 



l0g(l - Ml) 



log(l 



ui{ui 

U2) 



U2){ui 



- Us) 

log(l 



ys) 



dui du2 



U2{ui - U2){u2 - U3) ys{ui - ys){u2 - ys) 

For the map $, j = 3. The expression for T[log] appea ring in the stateme nt of the 
next theorem is complicated and is not given here, see |]Gil| , (1.1.31)] and |Gi2] for 
details. It depends only on the principal symbol (To(P) and involves certain Poisson 
brackets of the type {(Tq{B*), cro(P'')}, < t, s < 27r. 

For a function / e C°"{S*M) introduce the /cth Fourier coefficient along the 
closed geodesic of length 2n starting at a given point (x, ^) 



-ikt 



/(e*(x,0) 



dt 
2^' 



/c e Z. 



For an arbitrary B e 5'°(M), let us write 60 := cro(P), 6sub '■— sub(P), 6q(x,^) := 
&o(0*(a;,O)j s-iid omit the argument (x,^) G S*M for brevity. 
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Theorem 1.2. Let M he a Zoll manifold of dimension d G N. Let A be defined 
by (lOI). Assume that ai{A)(x,£,) = ai{A){x,-£,) for all (x,^) G T* M . Let B e 
^^{M) be such that <Jo{B) is strictly positive and a certain symbolic norm of L — B 
is sufficiently small. Then \ogB E ^"(M) and the following holds 
(1.8) 

Tr log PnBBa - Tr P„ (log B)Pn 



i / ^fc(log6o)Jlog6o)_;,da;d^ 



1 



+ (d-l) 



^fc(logfeo)fc ibsnh/bo)_^, 

k=l 

^ /'27r /'27r /'27T 7( t j ' 

E^^ i Jo Jo ^^'^*"''^^"^^"^^'^[i°g](&^,6o,6S)££i^ 



T[log](5o) Uxde 



Remark 1.2. The existence ol a full expansion of the type (1.8) has been proven 



in G01|. Explicit formulas for the coefficients of the first first two terms in (1 



(and also of the logn term which is a part of Tr P„ (log i?)P„) have been obtained 
in jGOlj |G01 . 



Remark 1.3. The formula (1.7) for $[log] has a structure similar to the coefficient 
in the third asymptotic term in a Szego type expansion for convolution operators 
obtained by R. Roccaforte in j^. A version of the Bohnenblust-Spitzer combina- 
torial theorem is used in the proof in Q. A second order Szeg o type expansion 
for convolution operators was established by H. Widom in |W5] with the help of 
the usual Hunt-Dyson combinatorial formula (3.4). A full asymptotic expansion 
for convolution operators has been found in | W€ ] . 



Remark 1.4. For d = I, the square bracket in ( |l.8| ) disappears. Also T vanishes, 
because all the Poisson brackets vanish in this case (for each of the two cotangent 
directions the angle does not change and (Jo{B) is homogeneous of degree in ^). 
Theorem 1.1 follows. 



Remark 1.5. As we have noticed above, Theorem and Theorem 1.2 give lower 
order corrections to (L6) for ip{u) = logu, as n ^ oo. In the next section we find 
explicit formulas for the coefficients in the asymptotic expansion of Tr log Pni3P„, 
as n — cxD. These expansions are fourth order for d = 1,2,3 and third order for 
d > 4, as n — !■ oo. 



Proof of Theorem l.i: an outline. We start by expanding ip{u) = \ogu in a power 
series about the point 1. After that it suffices to prove (1.8) for i!{u) — for all 
m e N. 



Let us recall the method of | G02| . Let tt^, fee N, be the projection on the fcth 
eigenspace of the operator A and set tt/c :— for fc < 0. Then P„ = X]I-=i ^fc for 
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71 e N, and we set P„ :— 0, n < 0. Recall that because M is a Zol l m anifold, 
for an arbitrary B e *0(7\.f)^ the operator B* , t e R, defined by ( [l~5|) is 27r- 
periodic. Therefore one can introduce the Fourier expansion B — X^jez where 
Bj G *°(M), j e Z, is defined by 



Bj = TTk+j B 



TTfe = 



fc=l 



1 

2^ 



For m £ N and ji , • • • , j„i G ^ introduce the notation 

M,n{j) ■■= min(0, ji, ji + j2, ■ ■ ■ , ji H h j™). 

The following commutation relation is of central importance here 

(1.9) BjP,,^P^+^B,, neN,jeZ. 

Using ( |1.S| ) one moves all the projectors to the left in the expression 

{^PnBP'ti) — ^ ^ PnB j^P'tiB j^Pyi ■ • • P^B j^P^ 

obtaining PnB"''Pn plus a correction. This implies for all n e N 
Tr(P„BP„r -Tr P„B"P„ 

= - ^ Tr ((P„ - PnPn+ji ■ ■ ■ Pn+ji + ---+j„,)Bj^ 

M„A~3) + 1 



■B, 



(1.10) 



H ^ Tr(7r„+fcPj, 

ii+---+jm=o fc=0 



Next, it is important that for any G G ^'"(Af), Af being a Zoll manifold, there exists 
a full asymptotic expansion for Tr(7rfcG'), as fc ^ oo, see Lemma 1.3 below. This 
result is due to Y. Colin de Verdiere |CdV |. The coefficients in that expansion are 
certain Guillemin-Wodzicki residues. Recall that for any compact closed manifold 
M of dimension d G N the Guillemin-Wodzicki residue of a pseudodifferential 
operator G G 4'™(Af) of order m G Z is defined by 



Res(G) 



a_d{G){x,C)dxd£,. 



S'M 



For an arbitrary G G ^"(A/) denote 

Ri{G) ■.= Rcs{A-'^+'G), 



/ = 0,1,2,.-. 



Lemma 1.3. Let M be a Zoll manifold of dimension d G N. Assume G G 5'°(Af). 
Then for any = 0, 1, 2, • • • , there exists Cn < oo such that 



N 



(1.11) 



Tr(^feG)-^fc'*-i-'Pz(G) 



1=0 



<CNk 



d-2-N 



fc G N. 



See jCdMI and [ p02| . Appendix] for the proof of ( PH ). 
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Now assume that for so me K < oo w e hav e Bk = for |fc| > K (this is not 
assumed in the fuU proof in Gi2 |). Then ( I.IC ) and Lemma [iT^ imply 
(1.12) 

Tr(P„BP„)'" -Tr P„B'"P„ 



E 



3l+---+3m=0 



d-1 



A/,„(j) 



o'o(-Bji • • • Bj^) dxdS^ 



S'M 



(Jo (-Bji • • • ) dxd^ 



S'M 

sub(Pji • • • Bj^J dxd£_ 



Now the symmetrization argument from |s|, pl| , G02, G03, 02 1 comes into 
play. Note that the domain of summation over ji, • • • ,jm in (1-12) is symmetric 
with respect to the permutations of the j's. Also the factors involving 

<Jo{Bn---BjJ^ao{B,,)-...-ao{BjJ 

are symmetric. Therefore one can write 



Jl+---+jm=0 



(To (Pji • • • Pj„ ) dxd£_ 



S*M 



E 



- J2 M,n {Jr)) I '^0 {B,, ■ ■ ■ ) dxd(„ 



where Sm is the set of all permutations of ji, • • • , jm and := [jrn ■ ■ ■ , jr ^)- 
After an application of the usual Hunt-Dyson combinatorial formula (see (3.4) in 
Section 3) 

E^u^.)- EE ""^"''"r''^'"^ 



res„ 



reSm k=l 



all summations but one in , , ■ „ become free and can be carried out. 

Now we also see the two difficulties with the computation of the n''"^ term. The 
first difliculty is in computing the minimum raised to power 2. Here one needs a 
generalization of the Hunt -Dyson formula, see Theorem 3J . The second difficulty 
is that the factor 



/ sub(Pjj • • • Pj,„) dxd^ 

J S'M 



is generally speaking not symmetric with respect to the permutations ji, ■ • • , jm, 
and the straightforward symmetrization as above fails. Indeed, recall that 



sub(P,, • • • P,,J = sub(P,J n ^o(%) 



(1.13) 



fc=i 
1 

2i 



p^fc 



^ {ao(P,J,(To(P,,)} n '^o^^)- 



l<k<l<m 



where {•, •} stands for the Poisson bracket. We see that the first sum in ( 1.13] ) is 
symmetric, whereas the second one is generally speaking not (not even after the 
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integration over S*M). However, each of the m{m — l)/2 terms in the second sum 
on the right-hand side of (1.13) possesses a partial symmetry. More precisely, we 
are allowed to permute the principal symbols that do not enter the Poisson bracket. 
This fact can be utilized in a modification of the symmetrization procedure. The 
corresponding contribution is lengthy, this is where the functional T arises. □ 

Remark 1.6. We would like to mention here that in the computation of T one really 
needs the information on the set of values {MmUr)} ^-es counted with multiplici- 
ties, that is the Bohnenblust-Spitzer theorem (see Section 3 for its statement), and 
not just a formula for X]res„, i^mUrW foi' P = 1; 2. 



2. Explicit asymptotic formulas for logdetP„i3P„, as n oo. 



Theorem 1.1 and 1.2 give an expression for Tr log P„BP„ = logdetP„i3P„ as a 
sum of TrP„(logP)P„ and the two lower order corrections, as n ^ oo. We would 
like to compute the coefficients in the asymptotic expansion of log det P„i?P„, as 
n — > oo. 

An auxiliary asymptotic expansion for 



(2.1) 



TrP„GP„ = ^Tr(7rfeG), 



fc=i 



need for dimension d 



where G E ^'''(M), is given in Proposition 2.1 below. A corresponding result 
for the coefficients of n'^, n'^~^ and logn for d > 2 can be found in |G02, after 
Lemma 0.2]. To prove Proposition 2.1, we sum over fc = I,-- - ,rt in (1.11), see 
|Gi2| for details. The subtle point is the constant coefficient in (2.1), which we 
1,2. The terms of all orders in (1.11) 



and also the possible 

rapidly decaying term, will contribute to it. Let us therefore for d = 1, 2 make an 
additional assumption 



(2.2) 

and set, for all A; G N, 
(2.3) efe 



\RiiG)\ < oo 



+ CX3 



.XG) ■.^Tv{^kG)-Y,k^'^~' Ri{G). 



If (2.2) holds, the series in ( p. 3D (and also in (2.4) below) is absolutely convergent, 
and also for any TV e N there exists cat < oo such that \ek\ < CNk~'^ , fc e N. Set 



(2.4) 



C{G) :-^e.(G). 



fe=i 



Note that in ( p.4| ) and also in the proposition below there only appears the series 
^j^qP;(G). However we need the absolute convergence (2.2) in the proof of the 
remainder estimate. 

Let 7 denote the Euler constant and ^ the Riemann zeta function. 

Proposition 2.1. Let M be a Zoll manifold of dimension d G N. Let P„ be as 
above and assume that G G \E'°(Af). For d — 1,2, assume in addition that (2.2) 
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holds, and let C{G) he defined by {2A). Then the following holds, as n ^ oo, 
(i) for d^l, 



TvPnGPn = n-Ro{G)+ logn • i?i(G) 

+ (c{G)+jRiiG) + f2ai)Ri{G))+^-[^R,{G)-R2{G))+o[^), 



(a) for d — 2, 



Tr P„GP„ = • ^ Ro{G) + " • i?o(G) + i?i(G)) + logn • i?2(G) 

OO ^ 

+ (g(G) + 7i?2(G) + ^ C(0 Ri+iiG)) + o(-' 



/=2 



(^iii^ for d > 3, 



Tr P„GP„ = n'' • i i?o(G) + n''-' ■ i?o(G) + ^ i?i(G)) 

+ n-^-' • i?o(G) + i i?i(G) + ^ i?2(G)) + logn • i?d(G) + ©(n'^- 



Remark 2.1. We see that TrP„(logB)P„ in Theorem yj and |L2| (G = logB) 
contributes to the leading asymptotic term of order n'', and also to all lower order 
terms of order , j = d — 1, • • • , 1, 0, —1, • • • , and to the logarithmic term logn, 
as n — !■ oo. In the classical SSLT the situation is much simpler: logi? is just the 
Toeplitz matrix of the operator of multiplication by log 6, and so Tr P„(log = 
(2n + l)(k^)o. 

Now we are ready to state the two corollaries. 

Corollary 2.2. Let B e ^'''(S^) have a strictly positive principal symbol and be 
such that a certain symbolic norm of I — B is sufficiently small. Assume also that 
(pT3) holds. Then the following holds, as n ^ oo. 



logdetP„BP„ = ci-n + cw • logn + cq + c_i • - + o(-^), 



where the coefficients are the sums of the corresponding coefficients from Theo- 
and Proposition \2.i\(i). 



rem 



1.1 



Assume further that ao{B) and sub(P) do not depend on the direction of ^, that 
is ao{B){x,^) — bo{x) and sub(P)(a;,^) = bsuh{x) , for (x,^) G S*S^ . Assume 
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also that 6-2 = 0. Then the following holds, as n 
(2.5) 

dx 



log det PnBPn 



27r 



log 60 (x) 



2tt 



log n ■ 2 



2tt 



&sub(a:) dx_ 
bo{x) 27r 

00 

^fc(log &o)fc(log foo)_ 



, fc=l 



1 



^fc(log feo)fc(fcsub/6o) 



C7(log S) + 7 i?i (log S) + ^ C(0 Ri (log S) 



, fc=i 



27r r 



+ 



^0(2;) V ^o(a;) 



2 n 



2^ 



0( ^ 



where C{logB) is given by (2.4). 

Remark 2.2. In some simple cases, for instance for 6sub(a;) — ±^60(2;), the left-hand 
side in (2.5) can be computed explicitly. The coefficients of n , lo gn, and ^ on the 
right in (2.5) in these cases are as expected, see also Remark 2.7. 

Corollary 2.3. Let M be aZoll manifold of dimension d > 2. Assume that P„ 
and A are as in Theorem Li. Let B G ^'''(M) have a strictly positive principal 
symbol and be such t hat the symbolic norm of L — B is sufficiently small. For d — 2, 
assume in addition (2^). Then the following holds, as n ^ 00, 



(2.6) logdetPnBPn^C^f^ ■n'' 



id) 



-C 



d-1-'- ^-d-2-- ^C'iog-logn + 0(n 
where the coefficients are the sums of the corresponding coefficients from Theo- 
and Proposition \2.\(ii) or (Hi). Lf one counts the logarithmic term, this 



expansion is fourth order for d = 2, 3 and third order for d> A. 

Remark 2.3. The coefficients C^'^', C^l\, and also c\^^, d g N, have been found in 
GOl| , |G02[ . 

Remark 2.4. The most interesting coefficient in (2^) is the constant one, since one 
can think of exp c'^^ as of a regularized determinant of B € (M) , see |G01 , 02|. 
The sum 

00 

-iRdilogB) + V Cmi+d-i{\ogB) 



1=2 



will for all d G N be a part of Cq'^^ For d = 1, Corollary gives a full expression 



for C^^^. For d 
lengthy. 



2, Corollary 2.2 gives a full expression for Cg ' , which is quite 



Remark 2.5. Let us compare the result of Corollary 2.2 with a generalization of 
SSLT to the case of B being an operator of multiplication by a function b{x) having 
discontinuities which is due to H. Widom and E. Basor. In this case log b(x) also has 

discontinuities, and so the series X^fegz 1^1 l(log^)feP diverges logarithmically. The 
following third order asymptotic formula holds for the operator of multiplication 
by a piecewise function b{x) 



(2.7) 



logAct PnBPn = ai ■ 71 + 02 ■ \ogn + 03 + 0(1), 



00, 



LOWER ORDER TERMS IN SZEGO THEOREMS ON ZOLL MANIFOLDS 



11 



where ai as in (2.5), the coefficient 02 has been computed by H. Widom in [W2|, 
and the cons tant term as has been found by E. Basor in Q. Note that the matrix 
of B in (2.7) is still Toeplitz, the logarithmic order of the subleading term being 
due to a slower decay of the Fourier coefficients of b{x). In our case the matrix of 
the operator B e ^'"(S^) is not Toeplitz (see Remark ^.7[ ), and the logarithm comes 
from the contribution of sub(i?). 

Remark 2.6. It would be interesting to find a compact formula for the constant 
term in (2^). We me ntio n that the constant 03 in (2/7) found in Q has a form 
similar to the one in (2^). It contains a "finite" term and and an infinite series 



of certain integrals multiplied by the values of the Riemann zeta function at the 
points 3, 5, • • • . Interestingly, an "invariant" form of that series has been found in 
W4 |. It is written as a single integral involving the function 



^{x) — logr(a:;). 



This gives the hope that a similar formula can be found for the constant (2.5). We 
have done some computations trying to find the constant term in (2.E), and the 
function ^'(x) has been appearing there. 



Remark 2.7. The matrix interpretation of Corollary 2.2 is as follows. Assume 



for simplicity that B g 4'°(§^) is as in the second part of Corollary 2.2, that is 
aa{B){x,^) = boix) and sub(B)(x,^) 6sub(a;) |^|~\ for all (x,^ G S*S^. Assume 
also that 5_2 = ^-3 = • • • = 0. Let Bq and i?sub be the operators of multiplication 
by bo and fegub, respectively. Let D be the linear operator in L'^{S^) such that 



De 



ikx 



1 ^ikx 



0. 



1^1 >1 

fc = 0. 



Note that this is not a differential, but rather a smoothing operator of order 
— 1. There is known a correspondence between the classical PsDO's on the cir- 
cle and their discrete counterparts, see Q for details. By that correspondence, 
the zeroth order PsDO B we started with, equals Bq + BsuhD. Introduce two 

Toeplitz matrices, Bq := {(6o)j_fc}j,/cGZ and Bsuh ■= {(&sub)j_fe}j,feez- Set also 



D := 



- i 1 1 



diag(- ,3,2, 
is Bq + Bsuh ■ D. Finally, set P„ 



. Then the matrix representation of Bq-\- B^^hD 

diag(- • • , 0, 1, • • • , 1, 0, • • • ) (2n + 1 ones). We 
see that Corollary 2.2 gives a fourth order asymptotics of the determinant of the 
truncated matrix P„ • [Bq + i?sub ■ D) ■ Pn- 



Now we can reformulate the question of ffiiding the constant term in (2.5) in 
purely matrix terms. Drop the hats and the dots for brevity. Let Ci be a Toeplitz 
matrix that corresponds to the operator of multiplication by fesub/^Oj and let the 
matrix D be as above. Clearly, the matrices Ci and D do not commute. Assume 
that the matrix log(/ — CiD) is well-defined. The question is to compute the 
constant coefficient in TrP„log(/ — CiD)Pn, or which is the same, the constant 
coefficient in 



(2.8) 



TrP„log(/-i?i/2Cii?i/2)P„ 



00. 



As we have noticed in Remark 2.1, this question is trivial for a Toeplitz matrix T 
in place of Di/^CiD^/^. 
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3. Generalized Hunt-Dyson combinatorial formula 

In this section we state the generahzed Hunt-Dyson formula (gHD) and remind 
the reader the Bohnenblust-Spitzer theorem (BSt). After that we briefly explain 



how the gHD is derived ft'om the BSt in |Gi3 . We refer to [Gi3, Section 4] for the 



first step of an independent proof of the gHD. This first step is a generalization of 
F. J. Dyson's idea on which the proof of the usual Hunt-Dyson formula in |^ is 



based. It is also explained in |Gi3, Section 4] how to reprove the BSt starting with 
the gHD. See [Gil, Chapter 2] for the details of the independent proof of the gHD. 

We state the result for the maximum, for the corresponding result for the min- 
imum one should replace the positive parts with the negative parts. For a e R, 
n e N denote 

(a)+ :=max(0,a), (a)!^ := ((0)+)". 

Fix any m € N and assume ai, . . . , a™ G M. Let Sm be the set of all permutations 
T of the numbers 1, . . . , to. For each t = ) G Sm denote 

Introduce the notation 



(3.1) MJar) 



max(0, a,-^ , a,-j + , • • • , + • • • + flr^ ), j = 1, • • • ,m, 
0, J - 0. 



Fix any ] = !,■■■ ,m. For arbitrary 

fci > 1, ■ • • , % > 1, ki + ■ ■ ■ + kj — m, 
we introduce the notation 

fci(a^) := H h a^^^ 

kj{ar) :— arfcj+...+fc^._^+i + ■ ■ • + (iTk^+...+kj_i+kj ■ 

Each of ki{ar), ^ = 1, • • • , j, is a sum of fc; permuted variables out of a^-j^, • • • , so 
that each of the permuted variables enters exactly one sum. Recall that (^^ " ^ ) :— 

^^"'^ I denotes a multinomial coefficient, here n,j £ N, h, ■ ■ ■ ,lj £ N U {0}, and 
+ - ■ ■ + = n. We are ready to state the generalized Hunt -Dyson formula (gHD). 

Theorem 3.1. For an arbitrary power n £ N, an arbitrary number of variables 

TO G N, and for arbitrary oi, • • • , a™ G M, the following holds 

(3.3) 

^ minM ^ / „ X (fci(a.))'j (fc,(a.))'^ 

2^ 2^ ji 2^ 2^ U^ ... i j ki k, 

fcl + --- + fo,-=m liH hi,- = n 



LOWER ORDER TERMS IN SZEGO THEOREMS ON ZOLL MANIFOLDS 



13 



Remark 3.1. In the case n — 1 we obtain j = 1, fci = m, h — 1, and (3.3) 
becomes the usual Hunt-Dyson combinatorial formula (HD) after (4.8)] 



(flr 



(3.4) ^ ' ^ 

= (to - 1)! (fli H h a,„)+. 

Recall now the statement of the Bohnenblust-Spitzer theorem (BSt) ||s[ The- 
orem 2.2]. It asserts that for any m G N and arbitrary ai, • • • , G R, the set 
{Mm{aT)}TeSm contains the same numbers with the same multiplicities as the set 
of sums of positive parts of the sums of Oi , • • • , Om, arranged according to the cyclic 
representations of all to! permutations. This becomes clear if we consider a simple 
example. Let us choose to = 3 and any ai,a2,a3 G R. The symmetric group S3 
consists of six permutations that can be written via the cyclic representations as 

^3 - {(123), (132), (12)(3), (13)(2), (23)(1), (1)(2)(3)} 

In this case the BSt states that the set 

(3.5) { max(0, Octi, flo-i +aa-2,aa-i + O-cri + acra)} ^^g,^ 

contains the same numbers with the same multiplicities as the set 
{(oi + 02 + 03)+, (ai + 03 + 02)+, 

(3.6) (ai + 02)+ + (03) + , (fli + 03)+ + (02) + , (02 + 03)+ + (ai) + , 

(ai)+ + (02)+ + (03)+}. 

Note that a certain maximum of zero and accumulating sums of the permuted 
variables does not need to equal the element of the set on the right-hand side 
corresponding to the cyclic representation of that permutation. The statement of 
the BSt is merely that the whole multisets are identical. 



Proof of Theorem 3.1: an outline. By the BSt, for any m S N and arbitrary real 



fli, • • • , Om, the counterparts of the sets (3.5) and (B^) contain the same numbers 



with the same multiplicities. Therefore for any function / : M+ — )■ C, the sum 



Sres li-^mi^a)) equals the sum of /'s values over the counterpart of (3.6) 



It turns out that for the multinomial function f{t) = t™ for any n e N, a further 



computation can be performed which leads to the gHD, see |Gi3, Section 2]. 

Without going into details let us just mention here that when we apply the 
multinomial formula in the sum over the analog of (^^), some terms have zero 
powers. Recall that we want all factors to be present in the right-hand side of the 
gHD (^^). However it turns out that the terms having at least one zero power can 
be summed together, and their sum gives exactly the sum of the nth power of the 
"previous" maximum M„i^i{acr) over Sm. 

At this step we use the principle of inclusion and exclusion and the Cauchy and 
Cayley identities from the theory of partitions. □ 

Remark 3.2. The steps of the derivation of the gHD from the BSt can be reversed. 
After that having started with the gHD, we can conclude that for any monomial 



/(t) its sum over (B.f;) equals the sum over (3.6). Due to an additional linearity, 
this actually holds for an arbitrary po lyno mial f(t). Now using the polynomial 
interpolation we arrive at the BSt, see [|Gi3| , Section 4]. 
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Remark 3.3. The independent proof of the gHD in |GIil, Chapter 2] proceeds by 
induction on the power n e N. The base of induction is the usual HD {n — 
1). In the proof of the inductive step, the key canceUation of the highest power 
n of the maximum after taking a sum over Sm foUows from a generahzation of 
F. J. Dyson's argument from the proof of the usual Hunt-Dyson formula in see 
Gi3, Section 4]. The proof of the inductive step is however quite technical. 



Remark 3.4. In |RS|, the authors rediscover a version of the usual Hunt-Dyson 
formula starting from the much more powerful BSt. 
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